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Abstract. This paper is part of a long term program to Cauchy-characteristic 
matching (CCM) codes as investigative tools in numerical relativity. The approach has 
two distinct features: (i) it dispenses with an outer boundary condition and replaces 
this with matching conditions at an interface between the Cauchy and characteristic 
regions, and (ii) by employing a compactified coordinate, it proves possible to generate 
global solutions. In this paper CCM is applied to an exact two-parameter family 
of cylindrically symmetric vacuum solutions possessing both gravitational degrees of 
freedom due to Piran, Saner and Katz. This requires a modification of the previously 
constructed CCM cylindrical code because, even after using Geroch decomposition to 
factor out the z-direction, the family is not asymptotically flat. The key equations in 
the characteristic regime turn out to be regular singular in nature. 
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1. Introduction 



This paper is part of a long term program to Cauchy-characteristic matching (CCM) 
codes as investigative tools in numerical tools in numerical relativity. For a review of 
CCM see the article of Winicour[|J. The approach has two distinct features: (i) it 
dispenses with an outer boundary condition and replaces this with matching conditions 
at an interface between the Cauchy and characteristic regions, and (ii) by employing a 
compactified coordinate, it proves possible to generate global solutions. This paper is 
continuation of work investigating CCM in systems with vacuum cylindrical symmetry 
(as a prototype system containing only one spatial degree of freedom). In |2| we 
developed the necessary machinery and in || we applied the approach to the exact 
solution of Weber and Wheeler [ffl which has one degree of polarization and so only 
possesses one gravitational degree of freedom. An agreement with the exact solution 
was found with a maximum error better then 1 part in 10 3 . In the same paper we 
also investigated the propagation of Gaussian wave packets possessing two gravitational 
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degrees of freedom, but this work did not involve a check against any exact solution. 
Such a test is likely to be a more rigorous one because, unlike the Weber- Wheeler case, 
it would involve passing derivative information across the interface. 

One of the major problems in numerical relativity is that there are very few exact 
solutions known which can be used to test numerical codes. In this paper CCM is 
applied to an exact two-parameter family of cylindrically symmetric vacuum solutions 
possessing both gravitational degrees of freedom due to Piran, Safier and Katz ||. This 
family is somewhat unphysical because the rotational degree of freedom diverges at 
future null infinity. Nonetheless, it can be used to test the CCM cylindrical code more 
rigorously because it involves passing derivative information across the interface. The 
previous CCM cylindrical code was constructed specifically for solutions which (after 
using Geroch decomposition to factor out the ^-direction) are asymptotically flat and 
therefore can not be used for processing the Piran et al family. However, by working 
directly in terms of the Geroch potential || it proves possible to develop a modified 
version of the code which is able to process the family. In section |2| we present the family 
of solutions and briefly discuss its properties. In section |3| we review the derivation of 
the field equations in the characteristic region for solutions which, after factoring out the 
z-direction, are asymptotically flat. In section ^ we discuss the asymptotic limit of the 
Piran et al family and derive the modified field equations and the interface equations. 
An appendix discusses the regular singular nature of two of the resulting equations. In 
section |5| we briefly discuss the results of the modified code. 

2. The two-parameter family of solutions 

Piran et al derive their metric by an indirect method starting from the Kerr line element 
in Boyer-Lindquist coordinates (t, r,9,<fi), namely 



ds 2 



A 



-(dt-a sin 2 8d<j)) 2 H — [(r 2 + a 2 )d<j) - adt} 2 + --dr 2 



+ P 2 d6 2 ,(l) 



P 



where 



2 ~2 i ~2 2 n 

p = r + a cos 0, 
A = r 2 — 2mr + a 2 . 

They first transform to a new isotropic radial coordinate R where 



(2) 
(3) 




(4) 



and then to cylindrical coordinates (t, p, Z, <p) where 



P 



Z = RcosO. 



RsinO 



(5) 
(6) 



Employing the complex trick of sending 



t — > iZ, p — >• p, Z — > it, (f) — > 0, a — >• ia, 



(7) 
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they then introduce new cylindrical coordinates (t, r, z, 0) given by 



t = - 

a 

P 



1 + 



m 2 + a 2 



4(p 



;,2 



t 2 



m 2 + a 2 



a [ 4(p 2 - t 2 ) 
z = Z — 2a~ 1 m{m + \J m 2 + a 2 )(f>, 

4> = QL(j), 



where 



a 



a/ m 2 + a 2 



(8) 

(9) 

(10) 
(11) 



(12) 



and the resulting line element is in the Jordan-Ehlers-Kompaneets form for a 
cylindrically symmetric spacetime M, 



3 27-2i/>/ 



(dt 2 - dr 2 ) + e 2 ^(dz + ud<j>f + e-'^r 



2 <M, 



If we finally introduce coordinates 

u = t — r, 
v = t + r, 

then we can express the two-parameter family of solutions in the form 



ds 2 

where explicitly 

e 2 ^ 

e 2 ^ 
u 

with 



-e 2 ^dudv + e 2lp (dz + ujdcf)) 2 + 

(X u + X v ) 2 + a 2 (l-X u X v ) 2 

(1 + A 2 )(1 + A 2 ) 
a 2 (l-X u X v ) 2 + (X u + X v ) 2 
a 2 E 2 + (X u - X v ) 2 



-2ip 



(v — u)' 



(13) 

(14) 
(15) 

(16) 

(17) 
(18) 



2a v / « 2 - 1 



a 



(a + Va 2 - 1) 



Wa 2 - TE(Xu + X l 



VX~X[a 2 (l-X u X v ) 2 + (X u + X v ) 2 ] 



X u 
X v 



1 + X U X V + 2a Va 2 - ly X U X V , 
(^/ a 2 + u 2 —uj /a. 



(19) 

(20) 
(21) 
(22) 



( V a 2 + v 2 + v) I a. 

The coordinates u and v are null, <fi is the canonical azimuthal coordinate and z lies along 
the axis of symmetry. The two parameters are a (0 < a < oo) which is a length scale 
and a (1 < a < oo) which is a measure of the total energy of the system, with a = 1 
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corresponding to flat space ||. The solution is regular everywhere in the coordinate 
range 

— oo < u < oo, (23) 

-oo < v < oo, (24) 

< < 2tt, (25) 

— oo < z < oo. (26) 

There is no conical singularity on the axis of symmetry and the solution reduces to 
Minkowski spacetime at past and future infinity However, the solution is conical at 
spatial infinity and singular at both past and future null infinity. More precisely, at 
future null infinity (cf Eq. (10) in ||) 

e*-l, ^ 7 - 1 T t r^' "--oo, (27) 

as v — > oo, with analogous behaviour at past null infinity on interchanging u and v. 
(Note that Eq. (10) in || includes a typographical error and that u and v should be 
transposed). 

Since Piran et al only derived this family of solutions by the indirect method 
described above, we tried to confirm that the family is vacuum by a direct computation of 
the Ricci tensor using both of the computer algebra systems SHEEP || and GRTENSOR 
0. Unfortunately, neither system was able to complete the calculation because the 
metric involves nested radicals and algebra systems are notorious for the difficulties such 
quantities present. However, we were able to confirm the indirect derivation described 
above. 

3. The standard field equations 

The standard treatment of the Jordan-Ehlers-Kompaneets cylindrically symmetric line 
element ( |i~3|) in the Cauchy region 

to < t < t f , (28) 
< r < 1, (29) 

is described in ||. In particular, it is shown that the independent set of dynamical 
equations for the variables if}, u, Lf and L are 

tt =~L, (30) 
r 

u, t = - 2e-^Lt (31) 

Lt,t = -e^l^r - -ru >rr - 2rtp :T uJ ir ), (32) 

Lt = VVV + rtr ~ ~e*y p + 2e-**(Lt)\ (33) 
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where Lf and L are defined in terms of the mixed components of the extrinsic curvature 
K»hy 

Lf = r 2 e^Kf, (34) 
L =r 2 e^{Kl-uKf). (35) 

These equations are augmented by the constraint equation 

7,r = - r^ r + \\l 2 + e-^(Ltn (36) 

which serves to determine 7 once the main variables are known. 

In the rest of this section we review the treatment in the characteristic region since 
this will need modification to cope with the Piran et al family of solutions. We introduce 
the compactified coordinate 

y = ~r = 1 » ( 37 ) 



in which case the line element becomes 

ds 2 = -e^-^du 2 + —e^-^dudy + e 2 ^(dz + uodcfy 2 + —-dtf, (38) 

v y 

where the metric functions ip, uj and 7 are now to be regarded as functions of u and y. 
The characteristic region then consists of 

u < u < Uf, (39) 
<y<l, (40) 

where the interface is taken to be at r = y = 1 and future null infinity is given by 
y = 0. Instead of working directly with the metric functions ip and u, we use the related 
quantities m and w which are defined by 

m = , (41) 

y 



o 

w — — 

y 



r p 1 f-p 

/ \ 2 {\y A uJ, y + yuj iU )du + - / \ 2 y 2 uj, y dy, (42) 
Jf y Jf 



where 



A = e 2 ^ = 1 + my, (43) 



o is the Geroch potential [0, and the integration is along any path connecting a fixed 
point F on the interface (we choose (u, y) = (ito,l)) to a general point P in the 
characteristic region (see Fig. [XJ) . With these definitions the vacuum equations can 
be written in the succinct form 

M =mj\, (44) 

W =wj\, (45) 

M 7 y = f(y, m, m t y, m m , w, w, y , W), (46) 

W tV = g(y, m, m >y , w, w tV , w m , M), (47) 
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Figure 1. Evaluation of the Geroch potential, o, at a general point P requires 
integration along the path F — > / — > P, see text for details. 



where / and g are explicit functions of their arguments [B|. Eqs. (44) and (45) serve to 



define M and W and ( |4"6| ) and fl4"T|) provide coupled propagation equations for M and 
W along the null rays ruling the hypersurfaces u = constant. The initial data consists 
of specifying m and w on the initial hypersurface {u = Uq, < y < 1} and the interface 
{uo < u < Uf, y = 1}, together with 7 at their intersection (u, y) = (uo, 1). The essence 
of the iterative integration scheme is to use fl46[) and ( ^7]) to determine M and W on the 



interior null hypersurface u = Ui and then fl44f) and ( ^5|) to determine m jU and w jU on 



u = Ui, which in turn determines m and w on the next neighbouring null slice. Finally 
7 is determined from the constraint equation (Eq. (42) in H) which has the form 

7,5, = /i(y, m, m >y , w, w iy ), (48) 

where h is an explicit function of its arguments. The CCM method requires the exchange 
of values of the metric functions and their derivatives at the interface at each iteration, 
the details of which are given in ||. 

4. The modified field equations 

The approach described in the last section breaks down for the Piran et al solution in 
the characteristic region because oj ^ — 00 as ?/ - >0. Let us consider the behaviour of 
the various functions on the hypersurface u = Ui in the asymptotic limit y — > 0. We 
have, 

K = (yfa 2 + u?-Ui)/a, (49) 

A„ = K/y 2 , (50) 



3/V, (51) 



where we set 



A„ = 2 + Ui y 2 + JA + Au lV 2 + (a 2 + u^y* /a = A/a + 0(y 2 ), (52) 



CCM for a family of cylindrical solutions 7 

S = X U X V + 2y^{l - a^)X u X v + y 2 = AXJa + 0(y). (53) 
Thus, asymptotically, 

e 2 ^ = l + 0(y), (54) 

= (l + a 2 Xl)/(l + Xl) + 0(y% (55) 



and 



(2J{a 2 -l)a\ u \ 1 2a v / ^^T/ ^ n ^, N , , 



\, J +0(y), (57) 
1 + cr Af, 



b 2J(a 2 - l)aX u 

™=- + V ' +0(y), (58) 
y 1 + a 2 A 2 



2J(a 2 - l)aX u 

b+ l + a 4 y+0{y) ' (59) 



where 



b = o{ui,0). (60) 

It is clear that A„, H, and tu are all divergent as y — ► 0. However, the ancillary 
quantities A„ and S as well as m and o, the Geroch potential, are all regular. This 
suggests rewriting the system (fH])-(l4"8|) in terms of these last two variables. We find 
explicitly that 

M = T ^- ) (61) 
1 + my 

O — j (62) 

1 + my 



M, = - ,°' v - O + ' 



y(l + my) 4(1 + my) 
1 + my 



-y(m + y 2 m t yy + 3ym )2/ ) 



+ {m 2 + 2ymm v + y 2 m 2 - o 2 



(63) 



° v =-° + TTT^ M ~ An V " ^ y°,w + °.») 
y (1 + my) 4(1 + my) 



+ 2(1 + my) 2 (m °' y + (64) 



7* = - 8 (i + mj/ ) 2 (^ + ^) 2 + ° 2 J- ( 65 ) 



Although the equations (|5"3| ) and ( |54] ) are coupled singular equations for M and O, 
they are regular singular equations and the solutions remain regular as y — > (see the 
Appendix). In fact, it is clear from the defining equations (Ell) and (|5"2"|), together with 
(|57|) and (|5T?D, that both M and O are of order unity in the limit y — > 0. These are 
the modified equations on which the new code is based for investigating the Piran et al 
solutions. 
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We need to augment the equations with the interface equations which relate the 
two sets of variables in the two coordinate systems on the interface. The relationships 
required for injection, i.e. going from the characteristic to the Cauchy region are, 

if) = -ln(l + rm/), (66) 
1 1 y 3 

i>,r = - ^yM - - — (ym),y, (67) 



^ rr = X -yM, a + l -y\yM), y - ^[( ym )^ 2 



1M 
27' 



+ ~{2(ym), y + [y(ym), y ] iy }, (68) 

(69) 



<*V = -?r, (70) 
y 2 \ 

O MO O l yO(ym) ty lyO, y 

U > rr - ~ y^X ~y\ X 2 P 2~T' 

L * = \v + \ y °> y - (72) 

Similarly, the relationships required for extraction, i.e. going from the Cauchy to the 
characteristic region are, 

A = e 2 ^, (73) 

m = r 1 ' 2 ^ - 1), (74) 



m 



y 



- r 2 e 2 ^(4V>,t + 4^, r + r- 1 ) + r, (75) 
m, w = 2r 7 / 2 e 2 ^(4^,« + 8^> + 8^ 2 + 10r~V,t + 16V>tVv + 4V> 

+ 8^ 2 + 10r _ V,r + r~ 2 ) - 2r 3/2 , (76) 

M = 2r~ 1/2 L, (77) 

o = f r'^Urdt, (78) 



F 











Ar^Lf - 2r 1 / 2 e^u r , (79) 
2r 2 e^(8ifj : tUJ ir + 2u; )tr + r _1 u; ir , + 8ip :r u :T + 2u :Tr + r -1 ^*) 

-8r 2 (Lt tt + Li >r ), (80) 

O = r~V^W (81) 



u 



</</ 



where the integration in (|78|) is along the interface from the point F : (t,r) = (t , 1) to 
the point / : (t, r) = (ti, 1) as shown in Fig. |l[ 

5. Results 

In order to test the CCM code against the Piran et al solution we will first compare the 
metric quantity if) in both the Cauchy and characteristic regions. Eq. (|lq) provides the 
exact solution in both cases when A n and A„, Eqs. ([H]) and (p2|) , are written in terms 
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of (t, r) coordinates and (u, y) coordinates for the Cauchy and characteristic regions 
respectively. It is not possible to use u for a similar comparison since from fl56|) we 
have seen that it diverges as y — > 0. We use, instead, the Geroch potential, o, which is 
computed directly by the code in the characteristic region, but must be constructed in 
the Cauchy region using its definition (cf Eq. (24) of 0), 

o = J r- x e^oo, r dt + j r~ l e^u )t dr. (82) 

It is straightforward to evaluate the integrals in the Cauchy region using finite-difference 
representations of u) >r and u >t . For the exact solution aj r and ou >t can be evaluated 
by differentiation, however the integral must be evaluated numerically to find a semi- 
analytic o. A similar approach is required for the 'exact' value of o in the characteristic 
region using Eq. (fEp. 

As mentioned previously, the two parameters a and a of the Piran et al solution 
represent the length scale and strength of the gravitational wave respectively. For t < 
the wave moves in from infinity, reaches its highest concentration at t — 0, where it 
rebounds from the r = axis, and is outgoing for t > 0. CCM allows us to evolve 
initial data containing an ingoing gravitational wave, something which is difficult to do 
for a Cauchy only evolution. Since we have set the interface at r = 1 the value of a can 
be chosen such that a substantial fraction of the the wave will move onto the Cauchy 
region and will be well resolved with the numbers of grid points we choose. Changing 
the interface position while keeping the 'interface distance to a' ratio the same results 
in an identical evolution on the Cauchy portion of the grid. 

The primary variable we used for comparison are ip and o (since 7 is a derived 
quantity). We considered grid resolution numbers of N = 301, 601 and 1201, where the 
error is defined as 

e(^) = ||^-^|| 2 /||^l|2, (83) 

where ipE is the exact value, tpc is the code computed value and || • H2 denotes the 
L 2 norm. We ran the modified code for a whole range of values in the a, a parameter 
space and found an error in ip of no more than 0.01% and an error in o of less than 
0.2%. However, in this version of the code, although the convergence rate of the solution 
starts as second-order it deteriorates to first order after long time evolution. In some 
recent work, colleagues in the Southampton Numerical Relativity group have modified 
the code so that it is now fully second order and very long time stable [|Kj. The 
improvement was achieved by dispensing with variables which are related by exponential 
or logarithmic functions at the interface, using the Geroch potential both in the Cauchy 
and characteristic regimes and also using an implicit method at the interface. This 
version of the code is currently being applied successfully to modelling dynamic cosmic 
strings. 

Surface plots are shown for evolutions of the metric quantity e 2 ^ in Fig. 
the Geroch potential o in Fig. [| and the radial derivative 7 >r (which indicates the 
distribution of energy within the wave)in Fig. f|, where the parameters are a = 0.5 and 
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Figure 2. Surface plot showing the time evolutions of the metric quantity e 2 ^ . The 
initial data has the parameters a — 0.5 and a = 10.0, representing a strong field 
spacetime. The interface is placed at z = 1 denoted by a gap in the surface plot. Note 
that for clarity the z-axis is reversed in this figure. 



a = 10. The surface plots use a radial coordinate z defined by, 

f r for < r < 1 „. 
z = { , ~ - (84) 

\2-yfor0<y<l V; 

so that there is a change in the coordinate system eA r — y — z — 1. Note that the value 
of e 2 ^ remains equal to 1 at null infinity, but is very close to zero at the axis r = 0. The 
incoming wave hits the axis at t = and rebounds. This can be seen very clearly in the 
surface plot of j jr . The large peak of Fig. 2c grows without bound as a is increased. 

For comparison purposes figures ^ |6] and shows the same quantities, but for 
parameter values a = 0.5 and a = 1.01, which represents an almost flat spacetime. 
Again the wave reaches its maximum concentration on the r = axis at t — 0, however 
the peak is much less dominant in this case. 
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Figure 3. Surface plot showing the time evolutions of the Geroch potential, o. 
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Appendix A. Regular-singular behaviour of modified propagation equations 

Let 



'41 



(1 + my) 



(A.1) 



which by ( 5?D and ([59]) is regular as y — > 0. The homogeneous part of the modified 
propagation equations (|63D and (|64] ) is 



F 

My+~0 

y 



o. 



O v -yFM - -O = 0. 

y 



(A.2) 
(A.3) 
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Figure 4. Surface plot showing the time evolutions of the energy distribution, 7,, 



Differentiating ( |A.3| ) with respect to y and using ( |A.3| ) to eliminate M and ( |A.2| ) to 
eliminate M tV we get 

O„-g + ^) O , + (* + ^ + ^) O = 0, (A.4) 
which has dominant singular part 

0,yy - -0, y + ^O = 0. (A.5) 

y y 2 

Substituting in the trial solution y k we obtain the auxiliary equation 

k 2 -3k + 2 = (A.6) 
which has roots k = 1, 2 and so gives rise to the regular independent solutions y and y 2 . 
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Figure 6. Same as Fig. |3|, but with parameters a = 0.5 and a = 1.01 representing an 
almost flat spacetime. 
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